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Optical trapping is a widely used technique, with many important applications in biology and metrology. 
Complete modelling of trapping requires calculation of optical forces, primarily a scattering problem, and non- 
optical forces. The T-matrix method is used to calculate forces acting on spheroidal and cylindrical particles. 
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1. Introduction 

Optical trapping provide three-dimensional 
confinement and manipulation of microscopic 
particles by a focussed laser beam. Optical trap- 
ping is a powerful and widespread technique, with 
the single-beam gradient trap (also known as op- 
tical tweezers) in use for a large number of bio- 
logical and other applications. 
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Figure 1. Schematic diagram of a typical optical 
tweezers setup 



The trapping beam applies optical forces (usu- 
ally divided into a gradient force, acting towards 
areas of higher irradiance, and scattering forces, 



including absorption and reflection forces) to the 
particle. 

The optical forces and torques result from the 
transfer of momentum and angular momentum 
from the trapping beam to the particle. Various 
approximate methods such as geometric optics or 
Rayleigh approximations are often used for the 
calculation of the optical forces. Such approx- 
imate methods are not necessary, since electro- 
magnetic scattering theory can be used for the 
calculation of forces, avoiding the limited ranges 
of applicability of the approximate methods. 

Other forces will also affect the motion of the 
particle. The most important of these forces, 
gravity, bouyancy, and viscous drag as the parti- 
cle moves through the surrounding fluid, are read- 
ily taken into account. 

2. Trapping as a scattering problem 

The optical forces and torques applied to the 
particle result from the transfer of momentum 
and angular momentum from the trapping beam 
to the particle. The total momentum transfer 
can be found by solving the electromagnetic scat- 
tering problem. A variety of numerical meth- 
ods can be used - finite element method, FDTD, 
discrete dipole approximation pQ, the T-matrix 
method |2l.ij . etc. 

A number of these have been used for opti- 
cal force calculations, including forces in optical 
traps | lloj . One method, however, stands out 
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as ideal for trapping force calculations - the T- 
matrix method. The T-matrix method can be 
considered an extension of Mie theory to arbitrar- 
ily shaped particles with arbitrary illumination. 
The main advantage of the T-matrix method 
is that trapping calculations usually involve re- 
peated calculation of the scattering for the same 
particle under differing illumination. In this case, 
the T-matrix need only be calculated once, since 
it is independent of the fields, whereas methods 
such as FEM, FDTD and DDA will require the 
entire calculation to be repeated. 

In the T-matrix method, the incident trapping 
field illuminating the particle is expressed as a 
sum of regular vector spherical wave functions 
(VSWFs): 
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jn(kr) are spherical Bessel functions, and dpm(^) 
are Wigner d functions. 

Similarly, the scattered fields are expressed as 
a VSWF expansion. In this case, since the far 
field must be an outgoing radiation field, 



( 2n + l 
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where M„ m (fcr) and N mn (fcr) are the same as 
RgM TO „(fcr) and RgN mn (fcr), with the spheri- 
cal Bessel functions replaced by spherical Hankel 
functions of the first kind, hn\kr). 

From the linearity of the Maxwell equations, 
there is a linear relationship between the incident 
and scattered fields: 

p = V T (11) a , ,+T (12) b , , (9) 
Fran / y J- mnm'n' "™ " ~ mnm'n' "m " \"J 
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The T-matrix can be calculated using the ex- 
tended boundary condition method (EBCM) 
For spherical particles, the T-matrix becomes 
diagonal, and the non-zero elements are the usual 
Mie coefficients. For rotationally symmetric par- 
ticles, the T-matrix is diagonal with respect to n. 
Computer codes to calculate T-matrices for such 
rotationally symmetric particles are available 6 . 

3. Representation of the trapping beam 

The use of the T-matrix method for scatter- 
ing calculations requires that the trapping beam 
be represented in terms of vector spherical wave 
functions, that is, the coefficients a mn and b mn 
in equation Q need to be found. The regular 
VSWFs RgM mn and RgN„ m provide a com- 
plete set of modes or partial waves, each indi- 
vidually satisfying the Maxwell equations, which 
can be used to represent any incident electromag- 
netic wave. For the simple case of a plane wave, 
E(r) = E exp(ik • r), with k in the direction 
(0,(f)), the expansion coefficients are |'2i:jj 

a mn = 4ir(-l) m i n d n C* mn • E exp(-im0) (11) 

b mn = 4f(-l) m i "- 1 <l„B t „ • E exp(-wn<£). (12) 

Note that the amplitude vector Eo contains the 
information regarding the polarisation and phase 
of the wave, and can be complex. 

In an optical trap, the incident field is usually 
a strongly focussed Gaussian or other beam. In 
principle, such a beam can either be decomposed 
directly into a VSWF representation, or into a 
plane wave spectrum, from which the VSWF ex- 
pansion coefficients can be found using equations 
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and i|12|l . In practice, this is problematic, 
since the usual descriptions of beams do not ac- 
tually satisfy the Maxwell equations. 

For the case of Gaussian beams, either plane 
wave expansion |7] or direct VSWF expansion 8, 
El can be used, although neither will give a beam 
identical to a traditional Gaussian beam. 

4. Optical forces 

Using the T-matrix method, with the T-matrix 
calculated by the publically available code by 
Mishchenko [H], and the beam shape coefficients 
in the localised approximation by Gouesbet |8I9| 
used to describe the beam, we calculated the vari- 
ation of the axial force acting on particles of vary- 
ing shape as a function of their position along the 
beam axis. 

The particles are polystyrene (n = 1.59) pro- 
late spheroids and cylinders, of varying aspect ra- 
tio as indicated (see figure |2J). The particles are 
of equal volume, with a volume equal to that of a 
sphere of radius 0.75/im, and are trapped in wa- 
ter by a Gaussian beam of waist width 0.8/im and 
free space wavelength 1064nm. 



of the beam, a negative force will act to axially 
trap the particle. If only optical forces are acting, 
the particle will come to rest at the zero optical 
force position where force curve crosses the zero 
force line with a negative gradient. 
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Figure 3. Axial force acting on spheroids of aspect 
ratios 1, 2, and 4. 




Figure 2. Differently shaped spheroidal and cylin- 
drical particles with aspect ratios of 1, 2, and 4. 



The trapping forces acting on the spheroid and 
cylinder with aspect ratio 1 is very small. This 
is due to interference due to reflections from the 
rear surface of the particle JT] . 

As the particles become smaller, the force be- 
comes less dependent on the shape of the par- 
ticle. Even for the particles considered here, it 
can be seen that the fine details of the shape (i.e. 
spheroid vs cylinder) only has a small effect on 
the force. As the particle becomes more elon- 
gated, it obstructs less of the trapping beam, and 
radiation pressure forces are reduced. 



The axial forces acting on the spheroids and 
cylinders are shown in figures |3| and 0] A nega- 
tive position on the beam axis indicates a position 
before the focal plane is reached, a positive po- 
sition is after the focus. A positive force acts to 
push the particle in the direction of propagation 



5. Non-optical forces 

A number of non-optical forces will act on the 
particle. Buoyancy and gravity are constant and 
are simply dealt with. Motion of the particle in 
the surrounding fluid is completely dominated by 
viscous drag. 



4 



3 
e 



s o 





a.r. = 2/ _\ 








a.r. = 1 












/ a.r. = 4 









Position on beam axis (nm) 



Figure 4. Axial force acting on cylinders of aspect 
ratios 1, 2, and 4. 



Since the time required for the particle to reach 
its terminal velocity of ~ 1/Ltms -1 is very short 
(r » 10~ 7 s is typical), we use r oc F instead 
of r oc F, since the particle will be moving at 
very close to the terminal velocity at all times. 
If the fluid is in motion, we use the velocity of 
the particle relative to the fluid. Since typical 
Reynolds numbers in trapping are extremely low 
(Re « 10~ 5 is typical), perfect laminar flow can 
be assumed - for a spherical particle, Stoke's Law 
will be an excellent approximation. 

In general, the trapping beam will heat the sur- 
rounding medium, possibly giving rise to convec- 
tive flow. While free convection is, in general, 
a difficult problem, the convection problem in 
trapping is perhaps the simplest possible. Due 
to the very small distances involved, the steady- 
state temperature distribution is reached in a very 
short time and a steady-state temperature 
distribution, independent of the convective flow, 
can be assumed. Similarly, any convective flow 
will reach steady-state very quickly. 

Further effects that can be included for a 
complete model include Brownian motion, ther- 
mophoretic effects due to uneven heating of the 
particle, the effects of nearby particles and sur- 
faces, etc. Accurate calculation of the optical 



forces can allow the evaluation of the accuracy 
of the modelling of these effects. 
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